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We study the processes pp — > W Z — > £v£ + £~ and pp — ¥ ZZ — > £ + £~i>v, where I = e or fj,. Using 
8.6 fb _1 of integrated luminosity collected by the DO experiment at the Fermilab Tevatron collider, 
we measure the WZ production cross section to be 4.50^Q'gg pb which is consistent with, but slightly 
above a prediction of the standard model. The ZZ cross section is measured to be 1.64 ± 0.46 pb, 
in agreement with a prediction of the standard model. Combination with an earlier analysis of the 
ZZ -> £ + ri + r channel yields a ZZ cross section of 1.44l^;|| pb. 

PACS numbers: 14.70.Fm, 14.70.Hp, 13.85.Q1 



I. INTRODUCTION 

This Article reports a study of WZ —> lv£ + £~ and 
ZZ — > £ + £~uv production in pp collisions at a center- 
of-mass energy of y 7 ! = 1.96 TeV. When not stated oth- 
erwise, we denote as Z the full Z/-f* propagator, with 
the requirement of 60 < Mu < 120 GeV on the dilepton 
invariant mass Mu for the decay Z — > £ + t~ . 

The pair production of W and Z gauge bosons probes 
the electroweak component of the standard model (SM) 
as these cross sections are predicted to be significantly 
larger in the presence of new resonances or anomalous 
triple-gauge-couplings [1, 2]. Diboson production is a 
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major source of background in many search channels for 
Higgs bosons. For example, ZZ decays correspond to 
some of the dominant backgrounds in searches for ZH 
production at the Tevatron. Understanding diboson pro- 
duction is therefore crucial for demonstrating sensitivity 
to the presence of a SM Higgs boson at the Tevatron. 

Production of W Z pairs was first reported by the CDF 
Collaboration, in 1.1 ftr 1 of integrated luminosity in 
the £v£ + £- channel 0. Evidence for WZ ^.re- 
production was also presented by the DO Collaboration 
in 1.0 ftr 1 of integrated luminosity 0]. The DO Col- 
laboration studied the same process with 4.1 ftr 1 [f|, 
measuring a production cross section of 3.90+J;gQ pb- 
The production of ZZ was first observed by DO in the 
£+£-£+£- final state with 2.7 ftr 1 [fj. Combination with 
an analysis of the £ + £~vv final state with 2.7 ftr 1 0], 
increased the significance from 5.3 to 5.7 standard devi- 
ations [f| . Evidence for ZZ production was also pre- 
sented by CDF in 1.9 ftr 1 8] of integrated luminos- 
ity. DO has recently analyzed 6.4 ftr 1 of integrated 
luminosity in the £ + £~£ + £~ final state Q. Combina- 
tion with the earlier 2.7 ftr 1 analysis [7] of the £ + £~vv 
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final state yielded a ZZ production cross section of 
1.40t°;^(stat)±0.14(syst) pb [icj|- The CDF Collabora- 
tion recently measured a cross section of 1. 64^38 pb f° r 
ZZ production using 6 fb^ 1 in the £ + £~£ + £~ and £ + £~vv 
final states [ll|. The ATLAS Collaboration has recently 
studied the WZ -> lvf+l~ and ZZ -> pro- 
cesses in collisions at = 7 TeV using 1.1 fb _1 of 
integrated luminosity 12l. [l3|. 

Following the approach of the previous DO analysis [7J 
of the ZZ — > £ + £~vv signal, we measure the ratios of sig- 
nal cross sections relative to the inclusive Z cross section. 
This has the advantage of cancelling the uncertainty on 
the luminosity, while other systematic uncertainties re- 
lated to lepton identification and trigger efficiencies, are 
also largely cancelled. 

All selection requirements and analysis techniques are 
optimized based on Monte Carlo (MC) simulation, or 
on data in signal-free control regions. The data are ex- 
amined in the region expected for signal only after all 
selection criteria are defined through simulation. 



II. 



APPARATUS 



The DO detector [TMl has a central-tracking system, 
consisting of a silicon microstrip tracker (SMT) and a 
central fiber tracker (CFT), both located within a 1.9 T 
superconducting solenoidal magnet, with designs opti- 
mized for tracking and vertexing at pseudorapidities [17j | 
\r)\ < 3 and \rj\ < 2.5, respectively. A liquid- argon and 
uranium calorimeter has a central section (CC) cover- 
ing |tj| up to ps 1.1 and two end calorimeters (EC) that 
extend coverage to \r/\ « 4.2, with all three housed in 
separate cryostats. The inter-cryostat (IC) region (1.1 
< 1 77 1 < 1.5) has little electromagnetic (EM) calorimetry, 
and reduced hadronic coverage. Additional scintillating 
tiles covering the region 1.1 < \r}\ < 1.4 provide improved 
energy resolution for hadronic jets. An outer muon sys- 
tem, covering \rj\ < 2, consists of a layer of wire chambers 
and scintillation trigger counters in front of 1.8 T toroidal 
magnets, followed by two similar layers after the toroids. 



III. DATA AND INITIAL EVENT SELECTION 

The data used in this analysis were collected with the 
DO detector at the Fermilab Tevatron pp collider at a 
center-of-mass energy of y/s = 1.96 TeV. An integrated 
luminosity of 8.6 fb _1 is available for analysis, following 
application of data quality requirements. 

Unlike the previous DO analyses of these channels, we 
do not restrict the offline event selection to events satis- 
fying specific trigger conditions. Rather, we analyse all 
recorded data in order to maximise the event yields. The 
analyzed events are mostly recorded by triggers that re- 
quire one or two electrons or muons with high transverse 
momentum, px- 



Since both signal processes involve the decay 
Z — > £ + £~ , a natural starting point is to select an inclu- 
sive sample of dilepton events. In addition to the e + e~ 
and dilepton channels, the ZZ analysis makes ex- 

tensive use of the e^ps^ channel for verifying modelling 
of dominant backgrounds (mostly WW —> £ + v£~v). In 
all channels we require that there are two oppositely 
charged leptons with an invariant mass Mu between 60 
and 120 GeV. The regions 40 < M u < 60 GeV and 
Mu > 120 GeV are used as control regions. The two lep- 
tons are required to originate from a common vertex that 
is located within ± 80 cm of the detector center along the 
z axis, defined by the beam direction. 

We define three different qualities for electrons and 
muons and refer to electrons or muons satisfying the cor- 
responding selection criteria, discussed below, as loose, 
medium, and tight, respectively. Electrons are treated 
differently when they are reconstructed in the CC, EC, 
and IC regions of the calorimeter. Loose CC/EC electron 
candidates are reconstructed using EM energy clusters in 
the calorimeter that satisfy minimal shower shape and 
isolation requirements and that are spatially matched to 
central tracks. A boosted decision tree (BDT) [l8[ is 
trained to separate electrons from jets, based on tracking, 
isolation, and shower shape information. Medium and 
tight CC/EC electrons must satisfy increasingly strin- 
gent requirements on the output from this BDT. 

In the IC region, there is no dedicated reconstruction 
of electrons. However, electrons traversing this region are 
identified using an algorithm that searches for hadronic 
decays of tau leptons. A neural network is trained to 
use shower shape, isolation, and tracking information 
to recover electrons from reconstructed taus while re- 
jecting hadronic jets. Electrons found in the IC region 
must be matched to a central track with at least one hit 
in the SMT and ten hits in the CFT. Loose, medium, 
and tight IC electrons must satisfy increasingly stringent 
requirements on the neural network output. In addi- 
tion, medium(tight) IC electrons must satisfy Ifrk/Pr 
< 0.2(0.1), where is the scalar pt sum of all tracks 
within an annulus of radius 0.05 < A1Z < 0.4 [l|| around 
the candidate electron. IC electrons are placed into two 
sub-categories: type-2(l) IC electrons do (not) contain a 
cluster of energy in the EM layers of the calorimeter. For 
type-2 IC electrons the momentum is determined from 
the calorimeter energy, whereas for type-1 IC electrons, 
we rely on the central track, which provides a relatively 
poor momentum resolution. 

Loose muons are required to have a track segment that 
has wire and scintillator hits in at least one layer of 
the muon spectrometer and a spatially matched track 
in the central detector. Loose muons must also sat- 
isfy a calorimeter isolation requirement of I™ on /pt < 
0.4, where 2^" on is the scalar sum of transverse ener- 
gies of all calorimeter cells within an annulus of radius 
0.1 < A7Z < 0.4 around the candidate muon. A track 
isolation requirement, l™° n /pT < 0.4, is also imposed 
on loose muons, where is the scalar px sum of 
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all tracks within a cone of radius A1Z < 0.5 relative to 
the candidate muon. Medium(tight) muons must satisfy 
l™f on /p T < 0.2(0.1) and l^ on /p T < 0.2(0.1). 

The number of events that pass an inclusive dilepton 
selection is used as the denominator for the purposes of 
measuring the ratio of WZ and ZZ cross sections to the 
Z cross section. The W Z dilepton selection requires two 
opposite charge medium quality leptons of the same fla- 
vor. Hard and soft electron(muon) px requirements are 
defined as p T > 20(15) GeV and p T > 15(10) GcV, re- 
spectively. IC electrons are considered only if they satisfy 
Pt > 20 GeV. The e + e~ and channels require that 

both leptons satisfy the appropriate soft pt requirement, 
and that at least one lepton satisfies the appropriate hard 
Pt requirement. 

The ZZ dilepton selection requires tight rather than 
medium leptons, and also includes the e ± /i =F control 
channel. The lepton px requirements are the same as 
in the WZ analysis. Since a reliable energy/momentum 
measurement is needed to suppress background from mis- 
reconstructed Z —> £ + £~ events, type-1 IC electrons are 
excluded. In the e ± /i^ channel, electrons are consid- 
ered only in the CC and EC regions. Figure Q] compares 
the Mu distribution in data and simulation (described 
in Section IIV|) after the ZZ dilepton selection, prior to 
any additional requirements. The data are well described 
by the simulation apart from the region of large Mu in 
the dielectron channel where the simulation yields more 
events than the data due to a mis-modelling of the reso- 
lution tails. This is shown not to have a significant effect 
on the analysis. 

The Z — > £ + £~ selections used as denominators in the 
cross section ratio measurements include some further 
requirements that are detailed in Sections IVl and IVlTI 



IV. PREDICTION FOR BACKGROUND AND 
SIGNAL 

Background yields are estimated using a com- 
bination of control data samples and MC simu- 
lation. The signal processes and certain back- 
grounds (WW_ £ + v£~v, ZZ ->■ £+£-£+£-, Z -> £+£~, 
tt £ + £~vvbb, Z~{ £ + £-~{ and Wj -> twy) are esti- 
mated based on simulations using the PYTHIA [2(| event 
generator and leading order CTEQ6L1 [2l[ parton dis- 
tribution functions (PDF). Events are passed through 
a GEANT ['T3 based simulation of the DO detector re- 
sponse. In addition, randomly triggered bunch crossings 
from data are added to the simulated events to model 
the effect of additional pp collisions. The GEANT based 
simulation does not include efficiency of the trigger. In- 
stead, the efficiencies of certain triggers (single-electron 
and single-muon) are measured using Z — > £ + £~ candi- 
date events from data. A parameterization of these effi- 
ciencies is applied to the simulated events, with a correc- 
tion determined from data for the additional efficiency 
that is gained from the remaining triggers (mostly dilep- 



ton and lepton-plus-jet triggers). The MC simulation 
does not accurately describe the dilepton p? distribution 
in Z production. Weights are therefore assigned to the 
simulated events as a function of their generated dilepton 
Pt, based on a comparison with more accurate predic- 
tions from RESBOS [23[ ■ The diboson events are similarly 
corrected as a function of diboson pt to match predic- 
tions from the next-to-leading order (NLO) event gener- 
ator powheg [13, Hl| . The simulation of WZ production 
by pythia does not include diagrams with 7* — > £ + £~. 
Weights are assigned to the generated events based on 
comparison of the Z — > £ + £~ invariant mass distribution 
with powheg, which does include these diagrams. The 
simulated events are further corrected for inaccuracies in 
reconstruction efficiency and energy/momentum resolu- 
tion for leptons. The MC predictions are normalized to 
match the observed event yield after the inclusive dilep- 
ton selection, such that knowledge of the integrated lu- 
minosity is not required. 

Instrumental backgrounds arise from the mis- 
reconstruction of hadronic jets as isolated electron and 
muon candidates. These backgrounds are estimated from 
data using the so-called "Matrix Method" [26j , since their 
rates cannot be modelled sufficiently well by the MC. We 
select a sample of events in which a high pt jet satisfies 
the trigger conditions for the event and is back-to-back 
in <j) with an electron or muon that satisfies the loose 
requirements. We measure the efficiency (ej c t) for a jet 
that is mis-reconstructed as a loose lepton to also satisfy 
the medium or tight lepton requirements. The equivalent 
efficiency for genuine electrons and muons (e s ig) is mea- 
sured in a sample of Z — > £ + £r candidate events. The 
Z+jets background in the W Z — > £v£ + £~ sample is esti- 
mated by selecting a sample of events in which the lep- 
ton assigned to the W — > £v decay is of loose rather 
than tight quality. Given the estimates of e s i g and £j e t, 
we solve a set of simultaneous equations to estimate the 
amount of background in the tight sample. A sample 
of Z+jets events generated with pythia is normalized 
to the estimate from data. The VF+jets background in 
the ZZ — ¥ £ + £~vv sample is estimated in a similar way, 
with a sample of MC events normalized to an estimate 
from data. In the inclusive dilepton sample, there is a 
small background from multijet events in the e + e~ chan- 
nel. This background is estimated by fitting the observed 
Mu distribution with the sum of simulated contributions 
and a multijet template that is obtained by inverting the 
electron quality cuts in real data. 



V. SELECTION OF WZ CANDIDATES 

Four channels are considered for WZ decay: e + e _ e ± , 
e + e^/i ± , n + fi^e^, fi^n - ^. Events must contain two 
oppositely charged medium quality leptons satisfying the 
Pt requirements described earlier and with Mu between 
60 and 120 GeV. The selection of WZ candidates fur- 
ther requires an additional electron (CC or EC) or muon 
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FIG. 1: Comparison of data and simulation in the Ma distribution after selecting an oppositely charged pair of tight quality 
leptons in the (a) e + e~ , (b) fi + pT , and (c) e*/^ channels. The lower halves of the plots show the ratio of data to simulation, 
with the yellow band representing the systematic uncertainty (see Section rVHI|) on the simulation. 



with pt > 15 GeV, and tight quality. This lepton must 
originate from a common vertex with the lepton pair 
that is assigned to the Z — > £ + £~ decay. If there are 
three like flavor leptons, there are two possible combina- 
tions of opposite charge leptons. In this case, the pair 
with smallest \Mu — mz\, where mz is the Z boson 
mass [27[, is assigned to the Z — > £ + £~ decay. Simu- 
lation shows that this assignment is correct in 93% of 
cases in the three electron channel, and 87% of cases in 
the three muon channel. In order to suppress the contri- 
bution from ZZ — > £ + £~£ + £~ decays, no additional re- 
constructed EM clusters are allowed for the £ + £~e ± se- 
lection, and no additional reconstructed muons for the 
£ + £~ /i^ selection. The additional EM clusters must sat- 
isfy Et > 5 GeV. Clusters that are not matched to a cen- 
tral track must satisfy loose shower shape requirements. 
There are no explicit px or isolation requirements im- 
posed on the additional muons. Events that satisfy these 
requirements, excluding the requirement of a third lep- 
ton, are considered as Z candidates, to be used in the 
denominator when measuring the ratio of WZ and Z 
cross sections. We choose to include a veto on more than 
two leptons in the Z selection, such that uncertainties 
in the veto efficiency are mostly cancelled in the ratio of 
WZ to Z cross sections. 

The WZ events are characterised by large missing 
transverse energy, ]fir , defined as the magnitude of a vec- 
tor sum of the transverse energies of cells in the calorime- 
ter. The IpT is corrected for muons, which only deposit a 
small fraction of their energy in the calorimeter, and for 
the energy loss of electrons. The variable Jftj, is defined 
by recalculating the $r after rescaling the transverse mo- 
menta of the leptons that are assigned to be Z daughters 
within 3 standard deviations of their resolution a(p^ ), 
through a fit that minimises the \ 2 function: 



where Tz is the total width of the Z boson [27J, and 

8pj) is the amount by which the pt of lepton i is shifted. 
This kinematic constraint is only used for the purposes 
of calculating the variable tfjj,- The requirement of Jf,^ 

> 20 GeV is imposed in order to reject Z+jets and Zj 
backgrounds. A background to the sub-channels with a 
W — > ev decay is the radiation of a high pr photon from 
a lepton in a Z — > £ + £~ decay. We therefore require that 
\Mm - 9 1.2 1 > | Mu - 91-2 1, where M m is the invariant 
mass of the three leptons. In addition, at least one of the 
leptons from the Z — ¥ £ + £~ decay is required to have px 

> 25 GeV. 

Tables U and [TT| list the observed and predicted event 
yields after all WZ selection requirements. The yields 
are also listed for the samples that exclusively fail the ]ftr 
or Mu requirements, but satisfy all other requirements. 
Figure [5] shows the Mu, and Mjf distributions for 
W Z — > £v£ + £~ candidate events, before imposing the re- 
quirement on each variable. The transverse mass is de- 
fined as = yj2pT^$r(X — cos Acf>), with pt being the 
transverse momentum of the charged lepton that is as- 
signed as the W daughter, and A<j> being the azimuthal 
angle between this lepton and the Jpr vector. Figure [3] 
shows the distributions of various kinematic quantities 
after combining the four sub-channels. 



VI. MISSING TRANSVERSE MOMENTUM 
ESTIMATORS 

The basic signature for the ZZ — > £ + £~vv analysis is 
a pair of charged leptons with invariant mass close to 
mz, produced in association with significant imbalance in 
transverse momentum, $t, due to the neutrinos from the 
Z — » vv decay. A substantial background corresponds 
to inclusive Z — > £ + £~ production in which either the 
leptons and/or any hadronic recoil is mis-reconstructed. 
Stringent selection requirements are necessary to sup- 
press this background, since (i) the production cross sec- 
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TABLE I: WZ — > £u£ + £ selection: Predicted and observed yields in the Z — > e + e sub-channels. The systematic uncertainties 
are provided for the predictions. 



e + e e + e ji? 

Accepted fail $r fail Mu Accepted fail Ipr fail Mu 



z ->■ i + r 


0.3 ± 0.1 


9 ± 1 








± 








3 ± 1 




7 


± 


2 







1 


± 








Z-y -> £+£"7 


0.6 ± 0.2 


10.1 ± 0.6 








± 








0.1 ± 0.0 





1 


± 














± 










0.6 ± 0.1 


1.0 ± 0.1 








± 








1.5 ± 0.0 





7 


± 











1 


± 








ti £+£~vvbb 


0.0 ± 0.0 


0.0 ± 0.0 








± 








0.0 ± 0.0 








± 














± 








ZZ -> i+Tvv 


0.0 ± 0.0 


0.0 ± 0.0 








± 








0.0 ± 0.0 








± 














± 








WW £+^"^ 


0.0 ± 0.0 


0.0 ± 0.0 








± 








0.0 ± 0.0 








± 














± 








Predicted background 


1.5 ± 0.4 


20 ± 1 





1 


± 








5 ± 1 




7 


± 


2 







2 


± 





1 


->■ ii^+r 


9.9 ± 0.2 


1.7 ± 0.1 





1 


± 





1 


13.9 ± 0.3 


2 


2 


± 





1 





6 


± 





1 


Predicted total 


11.4 ± 0.4 


21 ± 1 





2 


± 





1 


19 ± 1 




9 


± 


2 







8 


± 





2 


Observed 


17 


32 













17 






6 










1 







TABLE II: WZ — > lv£ + £ selection: Predicted and observed yields in the Z — > fi sub-channels. The systematic uncertain- 
ties are provided for the predictions. 





Accepted 


fail $ T 


fail M u 


Accepted 


fail $r 


fail M u 


Z -> £+£- 


1.5 ± 0.4 


12 ± 2 


0.5 


± 


0.2 


4 ± 2 


3 ± 1 


0.1 


± 


0.5 


Zj ->■ £ + r-y 


1.6 ± 0.4 


13.0 ± 0.5 


0.3 


± 


0.1 


0.1 ± 0.0 


0.1 ± 0.0 


0.0 


± 


0.0 


ZZ £+£-£+£- 


0.9 ± 0.2 


1.5 ± 0.2 


0.1 


± 


0.0 


1.6 ± 0.0 


0.7 ± 0.0 


0.1 


± 


0.0 


ti -> l+Tvvbb 


0.3 ± 0.0 


0.0 ± 0.0 


0.1 


± 


0.0 


0.1 ± 0.0 


0.0 ± 0.0 


0.0 


± 


0.0 


ZZ -J- i+t'vv 


0.0 ± 0.0 


0.0 ± 0.0 


0.0 


± 


0.0 


0.0 ± 0.0 


0.0 ± 0.0 


0.0 


± 


0.0 


WW £ + v£-v 


0.0 ± 0.0 


0.0 ± 0.0 


0.0 


± 


0.0 


0.0 ± 0.0 


0.0 ± 0.0 


0.0 


± 


0.0 


Predicted background 


4.3 ± 0.8 


26 ± 2 


1.0 


± 


0.3 


6 ± 2 


4 ± 1 


0.2 


± 


0.5 


WZ £u£+r 


14.0 ± 0.2 


2.1 ± 0.1 


0.9 


± 


0.1 


15.1 ± 0.4 


2.0 ± 0.1 


0.3 


± 


0.1 


Predicted total 


18.3 ± 0.8 


29 ± 2 


1.9 


± 


0.4 


21 ± 2 


6 ± 1 


0.5 


± 


0.6 


Observed 


26 


23 




3 




25 


12 
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tion for Z bosons exceeds that of the signal by four orders 
of magnitude and (ii) the rates for mis-reconstruction are 
difficult to simulate. Rather than estimate the genuine 
p'x hi the event, we use the approach of the previous 
DO analysis of this process [7], and construct variables 
that represent the minimum p'x consistent with the mea- 
surement uncertainties on the leptons and the hadronic 
recoil. 

First, the dilepton px is decomposed into two compo- 
nents with respect to a reference axis, t, as illustrated 
in Fig. We define t — p^ — p? 2 \ where p^ and 
Pj, ' are the px vectors of the two leptons. The dilepton 

Pt vector is defined as px = Px^ + Px^ > f rom which we 
define 



between the leptons, the a j; component is less sensitive 
to mis-measurement in the magnitude of the individual 
lepton transverse momenta than is off @, [H[ . For A0 < 
7r/2, this decomposition is no longer valid, and and 
aff are set equal to p 1 ^. 

The missing transverse momentum estimators, p" T , </l' Ti 
and fi L , are constructed as 



T 


= p^ - 




p T - 


h ^ ocoil - 


trkicts 

Vp T 


(5) 


T 


u 

= Wf - 


h2 


a s T - 


^recoil 


trkicts 

- a T 


(6) 


1 




h2 


4- 


^.recoil 
r a L 


trkicts 


(7) 



p 1 ^ — \pr\-, (2) The terms, p T , p^ coa and p^ k -' Gts (and similarly for ax 

HH _ | fi and cil) are corrections for lepton px mis- measurement, 

T ' ^ ' hadronic recoil measured in the calorimeter, and remain- 

a L = \Pt ' i\ > (4) ing hadronic recoil measured in the tracking system, re- 

spectively. These terms are described in more detail in 

where tisa unit vector in the direction of t. In the region the following sections. The factor of two is found to be 

A(j) > 7r/2, where A0 is the azimuthal opening angle optimal based on MC simulations. 
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FIG. 2: The distribution of (a-d) (e-h) Mu and (i-1) the W transverse mass of the WZ candidate events. The Jp^ 
requirement is not imposed for (a-d), and the Mu requirement is not imposed for (e-h). The rows correspond to different 
sub-channels as indicated on the figures. The vertical dashed lines indicate the requirements on ^ and Mu- The signal 
normalization is as described in Section Hvl 



A. Dilepton mis-measurement 

A correction for possible lepton px mis-measurement is 
determined by varying each individual lepton px within 
one standard deviation of its estimated uncertainty in 
order to separately minimise a^, off and p^ . Electrons 



that are reconstructed close to module boundaries in the 
CC or in the IC have relatively poor energy resolution 
and are given special treatment. The estimated uncer- 
tainty may be inflated to cover the difference between 
the calorimeter based pr measurement and the alterna- 
tive measurement from the central track. This is only 
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FIG. 3: Kinematic distributions for the WZ —¥ £v£ + £~ signal candidates after combining the different sub-channels. The 
following variables are shown: (a) the (b) the invariant mass of the Z — ► £ + £~ decay; (c) the W transverse mass; the 
transverse momenta of the (d) leading and (e) subleading leptons from the Z — > £ + £~ decay and (f) the charged lepton from 
the W decay; the transverse momenta of the reconstructed (g) Z ~ > £ + £~ and (h) W — > Iv decays. The vertical dashed lines 
indicate the requirements on ^ and Mu- The signal normalization is as described in Section [TV] 




Hadronic recoil 



FIG. 4: Illustration of the decomposition of the dilepton pr 
into or and cll components. 



allowed for the upward variation and protects against 
electrons for which the calorimeter has severely under- 
measured the energy. The amount by which, e.g., is 
reduced, is denoted a & T . These quantities are defined in 
such a way that they always carry a negative sign. 



Calorimeter recoil 



Two estimates of the calorimeter recoil are made, from 
the reconstructed jets and from the reconstructed $t- 
Jets are reconstructed using the DO mid-point cone al- 
gorithm with a cone size of AT?. = 0.5. They must 
be separated from the leptons by at least A7?. > 0.3 and 
satisfy px > 15 GeV. The px, &Ti ancl a L components 
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FIG. 5: Background vs. signal efficiency for (a) e + e and (b) fi + fJ. channels after varying the requirements on variables that 
are sensitive to the missing transverse momentum. The requirement <ji! T > 5 GeV is always applied. 



are calculated for each jet in the event, e.g., 

a j ; t(i) = P f (i, xi, (8) 

where pj, et ^ is the px vector of the ith jet. An individual 
jet that has a positive value (i.e., increases the momen- 
tum imbalance) is ignored. This approach ensures that 
jets which are not genuinely associated with the recoil 
system (e.g., from additional pp collisions or the under- 
lying event) are not allowed to generate a fake imbalance 
in an otherwise well reconstructed event. The sum of 
contributions from the jets is denoted, e.g., for the ax 

i jets 

component, a T . 

The Iffr estimate subtracts any contribution from the 
two leptons and then tests how well the remaining Ipx 
balances the dilepton system. Between the jet and ]pr 
based corrections we choose (separately for the ax, ah 
and pt components) the one that best balances the dilep- 
ton system. This correction term is denoted, e.g., a^ cml . 



C. Track recoil 

As a protection against events in which at least one 
hadronic jet fails to be reconstructed in the calorime- 
ter, we attempt to recover any remaining activity in the 
tracker. Track jets are reconstructed by merging together 
reconstructed tracks within cones of size A1Z = 0.5. 
These tracks must satisfy px > 1 GeV. Track jets must 
have at least two tracks within the cone, and be sepa- 
rated by at least A1Z = 0.3 from the leptons, and by at 
least A1Z = 0.5 from any calorimeter jets. Corrections to 
each of the (j>x, ax, and a£) components are determined 
in the same way as for calorimeter jets. 



D. Performance 

Figure [5] shows the Z — > £ + £~ background efficiency 
versus the ZZ —> £ + £~vv signal efficiency for a range of 
requirements on each of the variables; fi T , $ L , rf' T , and 
Ipx- The decays of Z — > t + t~ into e + e~, and 
e^/i^ final states produce a genuine missing px along 
the <zl direction. Our ZZ — > l + t~vv candidate selection 
requirements therefore include a "soft" requirement of </l' t 
> 5 GeV. The curves in Fig. [5] correspond to the combi- 
nation of this requirement and a varying requirement on 
the variable under study. The p" T variable has the best 
performance over the range of background efficiencies of 
interest. The efficiency for the requirement p" T > 30 GeV 
(and fi T > 5 GeV) is indicated explicitly by a star sym- 
bol. This is the requirement that is made in selecting 
ZZ — > l + l~vv candidates. The optimisation of the p" T 
requirement is discussed later. 



VII. SELECTION OF ZZ CANDIDATES 

Decays of ZZ into the following final states are con- 
sidered: e + e~vD, [l^h^vv. Events must contain two 
oppositely charged tight quality leptons satisfying the 
Px requirements described earlier and with Ma between 
60 and 120 GeV. To reject events in which the missing 
transverse momentum estimators defined in Section IVII 
are poorly reconstructed, we require that there are no 
more than two jets with px > 15 GeV and separated 
by at least A1Z = 0.3 from the leptons. In order to re- 
ject WZ -+ M+£- and ZZ -> £+£-£+£- events, there 
must be no additional EM clusters or muons according 
to the criteria in Section |V1 In addition, there must be 
no isolated tracks or hadronic taus with px > 5 GeV. 
This requirement is not necessary in the WZ analysis, 
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for which the ZZ — > £ + £~£ + £~ background is less signif- 
icant. These four types of objects are only considered if 
they are separated by at least A1Z = 0.3 from the lep- 
tons. The jet and additional lepton vetoes are also effec- 
tive in suppressing the background from tt — > £ + t~vvbb 
decays. The number of events that satisfy these require- 
ments is used as the denominator in the measurement of 
the ZZ/Z cross section ratio. Including the jet and lep- 
ton vetoes in the Z selection helps to reduce the impact 
of uncertainties in the corresponding efficiencies. 

Events are considered as ZZ — > £ + t~vv candidates if 
they further satisfy $ T > 5 GeV (to reject Z — > t+t~) 
and tf' T > 30 GeV (to reject Z —> e + e~ and Z —> 
Tables IIII1 IIVI and [V] show, for the three sub-channels, 
the predicted yields for each process. The yields are 
also presented for events that fail each requirement ex- 
clusively. Figure [6] shows the p" T and M« distributions 
before imposing their respective requirements. A neural 
network (NN) is trained to discriminate ZZ — > £ + £~vv 
from the dominant background in the final event sam- 
ple {WW — > £ + v£~v). The following input variables are 
used: the pr of each lepton, the ftp, the center of mass 
scattering angle cos 9* [3fJ] , the azimuthal angle between 
the leading lepton and the dilepton system Acj)(£i,££), 
and {Mm — mz)/cr{Mu) where cr(M«) is the estimated 
uncertainty on the measured dilepton invariant mass. 
Figure [6] also shows the NN output distribution of the se- 
lected signal candidate events. Separate NNs are trained 
for the e + e~ and channels, and the e + e _ ver- 

sion is shown for the e ± /Li =F channel. Figure [7] shows a 
number of kinematic distributions for the combination of 
ZZ — > e + e~vi> and ZZ — > fi + candidate events. 

Figure |8] shows how the predicted ZZ cross section 
measurement uncertainty varies as a function of the 
rf' T requirement. The expected systematic uncertainty 
rises rapidly below 25 GeV as the Z — > £ + £~ background 
starts to contaminate the sample. The requirement p" T 
> 30 GeV is close to the minimum and is in a region 
where the systematic uncertainty is small. 



VIII. SYSTEMATIC UNCERTAINTIES 

We measure the ratios of WZ and ZZ cross sections 
relative to the inclusive Z cross section. Lepton recon- 
struction, identification, and trigger efficiency uncertain- 
ties are largely cancelled in the ratio, as are those aris- 
ing from the vetoes on additional lepton candidates or 
other activity. The WZ analysis is sensitive to the lepton 
identification efficiencies, since the signal and Z — > £ + £~ 
samples differ by the requirement of an additional tight 
quality reconstructed electron or muon. The ZZ analysis 
is sensitive to the modelling of the diboson px, since re- 
quirements on the missing px estimators are less efficient 
in signal events with a large hadronic recoil. Tables IVII 
and IVIII list the sources of systematic uncertainty on the 
WZ and ZZ cross section measurements, respectively. 
We list the fractional variations in the number of pre- 



dicted background events iVb g d, the acceptances (mul- 
tiplied by efficiencies) for signal (A s ; g ) and Z — > £ + £~ 
{Aa), and the measured signal cross section er s i g . The 
following sources of systematic uncertainty are consid- 
ered. 

• Beam conditions 

The differential distributions of the instantaneous 
luminosity and vertex z position are varied to cover 
any disagreement with the data. 

• Physics modelling 

The value of the gi parameter in resbos is var- 
ied when determining the corrections that are ap- 
plied to the simulated Z — > £ + £~ events. This is a 
model parameter that describes the intrinsic trans- 
verse momentum of the partons within the collid- 
ing hadrons. As a test of sensitivity to the diboson 
Pt modelling, the reweighting in this variable is 
switched off. 

• Jet reconstruction 

The jet energy scale, resolution, and reconstruc- 
tion efficiencies are varied within their uncertain- 
ties. The simulation requires additional corrections 
to the energy response for jets in the IC region. 
An additional systematic uncertainty is assigned to 
these corrections. The track jet reconstruction effi- 
ciency is also varied to cover an observed disagree- 
ment with the data. 

• Lepton momentum scale and resolution 

The lepton momentum scales and resolutions are 
varied within their uncertainties, as are the re- 
construction and identification efficiencies. Non- 
Gaussian tails in the lepton momentum resolution 
are also considered. 

• Instrumental backgrounds 

The VF+jets and Z+jets background normaliza- 
tions are varied within the uncertainties of the esti- 
mate from data. All other variations on the simula- 
tion (e.g., lepton momentum scales and resolutions) 
are allowed to vary the shape of these backgrounds. 
Since PYTHIA does not include the matrix element 
for wide angle photon emission in W~f production, 
the normalization of this process is varied by a fac- 
tor of two, which is considered to be an overesti- 
mate but introduces no significant uncertainty on 
the ZZ cross section measurement. 

• Trigger efficiencies 

The trigger efficiencies are estimated to introduce 
a negligible uncertainty into the cross section mea- 
surements. 
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FIG. 6: (a-c) The rf' T distribution of the ZZ —¥ t + t~vu candidate events before imposing the rf' T requirement, (d-f) The Mu 
distribution of the ZZ — > i + i~vv candidate events before imposing the Mu requirement, (g-i) The neural network output 
distribution of the accepted ZZ — > l + l~uv candidate events. For the e fi^ channel, the neural network trained in the e + e~ 
channel is shown. The vertical dashed lines indicate the requirements on rf' T and Mu- The signal normalization is as described 
in Section ITVl 



TABLE III: Table of predicted signal and background yields for the ZZ — >• e + e vv signal and control regions. The systematic 
uncertainties are provided for the predictions. 

Rejected by requirement on 



Process 


Accepted 


fr 






4't 






Mu 


Extra lep. 


Charge 




Jets 




Z -> e+e" 


0.6 ± 0.3 


11666 ± 1665 







± 


1 




0.3 ± 0.2 


3 ± 2 


0.0 


± 0.0 





1 


± 





1 


Z T + T~ 


0.1 ± 0.1 


8 ± 2 


1 


4 


± 





2 


0.0 ± 0.0 


0.0 ± 0.0 


0.0 


± 0.0 








± 








WW -> i + vl~v 


35 ± 1 


35 ± 1 


1 


7 


± 





1 


33 ± 1 


9 ± 5 


0.3 


± 0.1 





1 


± 





1 


wz ->■ ivi+r 


2.3 ± 0.1 


1.9 ± 0.1 





2 


± 








0.2 ± 0.1 


14 ± 2 


0.2 


± 0.1 








± 








W — >• ev 


6 ± 2 


13 ± 2 





3 


± 





1 


5 ± 1 


2 ± 1 


4 


± 1 








± 










3.3 ± 0.3 


5.5 ± 0.5 








± 





1 


2.8 ± 0.5 


0.6 ± 0.5 


3.3 


± 0.4 








± 








zz -»■ i + £-£ + e~ 


0.0 ± 0.0 


0.1 ± 0.0 








± 








0.0 ± 0.0 


1.3 ± 0.2 


0.0 


± 0.0 








± 








tt -t l + l~wbb 


1.0 ± 0.2 


1.4 ± 0.2 





4 


± 





1 


1.2 ± 0.1 


7 ± 1 


0.0 


± 0.0 





2 


± 





1 


Predicted background 


48 ± 2 


11749 ± 1668 




4 


± 


1 




43 ± 2 


37 ± 11 


8 


± 1 





4 


± 





2 


ZZ -> l + l~vv 


13.6 ± 0.4 


7.4 ± 0.2 


1 


3 


± 





1 


0.6 ± 0.0 


4 ± 2 


0.2 


± 0.0 





1 


± 








Predicted total 


62 ± 3 


11756 ± 1668 




6 


± 


1 




43 ± 2 


41 ± 13 


8 


± 1 





4 


± 





2 


Observed 


61 


10560 






12 






50 


63 




12 






1 
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FIG. 7: Distributions of (a) cos(0*), (b) (Mm — mz)fo'(Mu), (c) IpT, transverse momenta of the (d) leading and (e) sub- 
leading lepton, (f) the azimuthal angle between the leading lepton and the dilepton system A<f>(£i,££) for the combination of 
ZZ — s> e + e~vu and ZZ — > fiTvv candidates. The signal normalization is as described in Section HVl 




Pj requirement (GeV) p| requirement (GeV) 

FIG. 8: Variation of the predicted uncertainties on the measured ZZ cross section with the choice of rf' T requirement in the 
(a) e + e~ and (b) fi + yT channels. 



IX. MEASUREMENT OF CROSS SECTIONS 

The ratios of the signal (WZ or ZZ) cross sections to 
the inclusive Z cross section are determined as follows: 

_ N^/(A sig x B sig x £) 
N$»/(Au x B u x C) ' 

where C is the integrated luminosity; Bu and B s - lg are the 
known branching fractions for Z £ + £~ and the signal 



decay, respectively [27]. We choose an acceptance win- 
dow of 60 < M u < 120 GeV. 

The number of observed signal events, N°^ s is deter- 
mined by allowing the predicted signal yield to float such 
that the following likelihood function is maximized: 

bins 

L = l[V(N° hs ;N? Ied ), (9) 

i=0 

where V is the Poisson probability to observe N° hs events 
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TABLE IV: Table of predicted signal and background yields for the ZZ —¥ vv signal and control regions. The systematic 

uncertainties are provided for the predictions. 

Rejected by requirement on 



Process 


Accepted 














Ma 


Extra lep. 


Char 


ge 




Jets 




Z — y fi + fi~ 


0.3 ± 0.5 


8519 ± 1372 




3 


± 


6 




2 ± 2 


3 ± 2 


0.4 


± 


0. 


1 


0. 


.1 


± 





4 




0.0 ± 0.0 


5 ± 3 


1 


.4 


± 


0. 


3 


0.0 ± 0.0 


0.1 ± 0.0 


0.0 


± 


0. 





0. 


.0 


± 





.0 


VV VV — 7 £ VX, V 


qi j. o 
ol zt Z 


A 8 -1- 9 
It Z 


1 
i 


.4 




u. 


o 
Z 


zy it z 


Q 4- 

y it o 


0.0 


± 


0. 





u. 


.1 


± 





1 

. 1 


wz -> ivi+r 


2.0 ± 0.1 


2.9 ± 0.2 





.2 


± 


0. 





0.3 ± 0.0 


12 ± 2 


0.2 


± 


0. 


1 





.0 


± 





.0 


W -> fxv 


2.3 ± 0.4 


9 ± 2 





.0 


± 


0. 


1 


2.9 ± 0.7 


1.1 ± 0.9 


0.8 


± 


0. 


2 





.0 


± 





.0 


zz -> e+re+r 


0.0 ± 0.0 


0.2 ± 0.0 





.0 


± 


0. 





0.0 ± 0.0 


0.9 ± 0.2 


0.0 


± 


0. 








.0 


± 





.0 


ti ->• t+Tvvbb 


0.8 ± 0.1 


2.0 ± 0.2 





.4 


± 


0. 


1 


0.9 ± 0.1 


6 ± 1 


0.0 


± 


0. 








.2 


± 





.1 


Predicted background 


36 ± 2 


8602 ± 1374 




6 


± 


6 




36 ± 2 


32 ± 9 


1.4 


± 


0. 


2 


0. 


4 


± 





.3 


ZZ -> l+l-vv 


11.8 ± 0.3 


11.0 ± 0.3 





.9 


± 


0. 


1 


0.8 ± 0.1 


4 ± 2 


0.0 


± 


0. 








.0 


± 





.0 


Predicted total 


48 ± 2 


8613 ± 1374 




7 


± 


6 




36 ± 2 


35 ± 11 


1.4 


± 


0. 


2 





4 


± 





.3 


Observed 


58 


7416 






7 






42 


60 




4 










1 







TABLE V: Table of predicted yields in the e /i+ channel. The systematic uncertainties are provided for the predictions. 

Rejected by requirement on 



Process 


Accepted 








4' T 






M u 


Extra lep. 


Char 


ge 




Jets 




Z -> e+e" 


0.0 ± 0.0 


17 ± 7 


0. 


.0 


± 





.0 


0.0 ± 0.0 


0.0 ± 0.0 





.0 


± 


0. 


.0 





.0 


± 


0. 


.0 


Z — > 


0.0 ± 0.0 


6 ± 2 


0. 


.0 


± 





.0 


0.3 ± 0.3 


0.3 ± 0.1 





.0 


± 


0. 


.0 





.0 


± 


0. 


.0 


Z T+r- 


0.1 ± 0.1 


19 ± 14 


4 


.5 


± 





.4 


0.1 ± 0.1 


0.2 ± 0.2 





.0 


± 





.0 





.0 


± 





.0 


WW -)■ £+v£~v 


69 ± 3 


84 ± 4 


3. 


.7 


± 





.2 


67 ± 3 


19 ± 11 





.4 


± 


0. 


.1 





.3 


± 


0. 


.2 


WZ -> lvl + l~ 


0.4 ± 0.1 


0.7 ± 0.1 


0. 


.0 


± 





.0 


0.4 ± 0.1 


3.3 ± 0.6 





.4 


± 


0. 


.1 





.0 


± 


0. 


.0 


W -> ev 


4 ± 2 


9 ± 3 


0. 


.1 


± 





.2 


5 ± 1 


1 ± 1 


1 


.2 


± 


0. 


.6 





.0 


± 


0. 


.0 


W -> flU 


5 ± 4 


12 ± 9 


0. 


.2 


± 





4 


5 ± 4 


1 ± 2 




3 


± 


2 







.0 


± 


0. 


.0 


W"f — > eu-y 


3.4 ± 0.5 


6.4 ± 0.4 


0. 


.1 


± 





4 


3.4 ± 0.3 


0.7 ± 0.6 


3 


.2 


± 


0. 


.3 





.0 


± 


0. 


.0 


ZZ -> i+ri+r 


0.0 ± 0.0 


0.0 ± 0.0 


0. 


.0 


± 





.0 


0.0 ± 0.0 


0.2 ± 0.0 





.0 


± 


0. 


.0 





.0 


± 





.0 


ti i+Tvvbb 


2.3 ± 0.2 


3.3 ± 0.2 


1. 


.0 


± 





4 


2.1 ± 0.2 


13 ± 3 





.0 


± 


0. 


.0 





.3 


± 


0. 


.1 


Predicted background 


84 ± 6 


157 ± 19 


9. 


.6 


± 





.5 


83 ± 6 


39 ± 16 




8 


± 


3 







.5 


± 


0. 


.3 


ZZ -> t+rvv 


0.0 ± 0.0 


0.0 ± 0.0 


0. 


.0 


± 





.0 


0.1 ± 0.0 


0.0 ± 0.0 





.0 


± 


0. 


.0 





.0 


± 


0. 


.0 


Predicted total 


84 ± 6 


157 ± 19 


9. 


.6 


± 





.5 


83 ± 6 


39 ± 16 




8 


± 


3 







.5 


± 


0. 


.3 


Observed 


73 


162 






7 






96 


60 






8 


















in the ith bin, given a prediction of iV^ 10 . In the WZ 
analysis, the Mt distribution is used, while the neural 
network output distribution is used in the ZZ analysis. 
The 68% C.L. interval on the signal yield is defined by 
(5(lnL) = 0.5, with respect to the maximum of InL. 

Table rvTTTI lists, for the six different sub-channels, the 
ratios of inclusive Z and signal acceptances that are es- 
timated from the simulation. Table IIXI lists the mea- 
sured 71 values. The p- values for consistency of the dif- 
ferent sub-channels are 54% and 11% for the WZ and 
ZZ analyses, respectively, evaluated using a \ 2 test. For 
the combination of respective sub-channels, we measure: 

K(WZ) = 0.593 ± 0.080(stat) ± 0.017(syst)(xl0" 3 ), 

K(ZZ) = 0.216 ± 0.058(stat) ± 0.017(syst)(xl0" 3 ). 

A theoretical calculation of the Z cross section can be 
used to translate these into signal cross section measure- 
ments. The product of the cross section and branching 



fraction for Z — > £ + l (one lepton flavor) is calculated us- 
ing a modified version of the next-to-NLO (NNLO) code 
of Ref. [1H with the MRST2004 NNLO PDFs [32[. Since 
this code excludes the 7* and Z/j* interference, a cor- 
rection factor is determined using PYTHIA and the NLO 
event generator MC@NLO [13]. For 60 < Mu < 120 GeV, 
the result is, 

a(pp -» Z/7*) x Bu = 255.8 t^o pb, 

where the uncertainties arise from variations in the PDFs 
and the renormalization and factorization scales, and 
with B u = 3.3658 ± 0.0023 % [13 . The measured WZ 
cross section with 60 < M« < 120 GeV is 

a{pp ->■ WZ) = 4.50 ± 0.61(stat)+°;25(syst) pb. 

This result is slightly higher than, but still consistent 
with a prediction of 3.21±0.19 pb from the NLO program 
MCFM [13 with the MSTW2008 NLO PDFs ||, and 
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TABLE VI: Table of uncertainty sources in the W Z cross 
section measurement after combining the four sub-channels. 
All values are given in percent. 



iinst profile 
Vertex z profile 
Z/-y*p T 
Diboson pt 
Jet energy scale 
Jet energy resol. 
IC jet treatment 
Electron pr scale 
Electron pr resol. 
Electron pr tails 
Muon pr scale 
Muon pr resol. 
Muon pt tails 
Track reconstr. 
Muon reconstr. 
Electron reconstr. 
Z/7*+jets model. 

Systematic 
Statistical 
Stat. syst. 



4.0 
1.6 
0.0 
0.1 
6.0 
2.2 
1.1 
0.3 
1.0 
0.1 
0.1 
0.9 
1.0 
0.1 
0.2 
0.2 
17.7 



3.3 
0.9 
0.0 
0.4 
0.3 
0.0 
0.0 
0.1 
0.0 
0.3 
0.1 
0.1 
0.4 
1.1 
0.5 
0.2 
0.0 



19.4 2.9 3.7 



19.4 2.9 3.7 



0.9 
0.4 
0.0 
0.4 
0.2 
0.0 
0.0 
0.1 
0.0 
0.4 
0.1 
0.0 
0.2 
0.3 
0.2 
0.0 
0.0 
1.2 



1.2 



0.2 
0.7 
0.2 
0.2 
1.3 
0.2 
0.2 
0.2 
0.2 
0.2 
0.2 
0.2 
0.2 
0.7 
0.2 
0.2 
2.5 
3.1 
13.2 
13.6 



setting the renormalization and factorization scales equal 
to mw + iTiz- The measured ZZ cross section with 60 < 
M u < 120 GeV is 



<j{pp —> ZZ) 



1.64±0.44(stat)+°;^(syst) pb. 



This can be compared to a prediction of 1.30 ±0.10 pb 
from MCFM setting the renormalization and factoriza- 
tion scales equal to 2m z- For comparing to and com- 
bining with previous measurements it is more convenient 
to correct the cross sections for the presence of diagrams 
involving 7*. A correction of +3.4% is obtained by com- 
paring ZZ cross sections with and without 7* and Z/j* 



interference from MCFM [34|. Combining this corrected 



ZZ cross section with a previous DO measurement Q in 
the ZZ -> £+£-£+£- channel yields 



cr(pp — > ZZ) 



= l.U±° H (stat) 



+0.17 
-0.19 



(syst) pb. 



X. 



CONCLUSIONS 



We measure the production cross sections for the pro- 
cesses pp — > WZ — > and pp — > ZZ — > vvl + l~ , 
using 8.6 fb _1 of integrated luminosity collected by the 
DO experiment at the Fermilab Tevatron collider. For de- 
cay channels involving electrons and muons, we observe 
agreement between the different sub-channels as can be 
seen in Fig. |H1 Combining the sub-channels yields a WZ 
cross section of 4.50 io.66 P^i which is slightly above, 



TABLE VII: Table of uncertainty sources in the ZZ cross sec- 
tion measurement after combining the e + e~ and chan- 
nels. All values are given in percent. 
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A 
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L inst profile 
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2 


0.7 


1.8 


Vertex z profile 


1 
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0.0 


0.6 


Diboson pt 


2 


6 








1 
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1.8 


3.7 


Jet energy scale 


1 


1 





8 


1 


5 


0.8 


1.8 


Jet energy resol. 





9 
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1 


0.0 


1.8 


IC jet treatment 





2 
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4 


0.2 


0.6 


Jet reconstr. 
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0.2 


0.0 


Trkjet reconst. 
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1.2 


3.1 


Electron pt scale 





4 














0.0 


0.6 


Electron pr resol. 


1 








1 





5 


0.4 


1.8 


Electron pt tails 


1 














6 


0.6 


1.2 


Muon pt scale 





1 














0.0 


0.0 


Muon pt resol. 





5 





1 





5 


0.5 


0.6 


Muon pt tails 





1 





1 





5 


0.4 


0.6 


Lepton eff. vs pt 




















0.0 


0.6 


Lepton eff. vs 77 




















0.0 


0.6 


VF+jets model. 


1 


9 














0.0 


0.6 


W'y model. 


3 


9 














0.0 


1.8 


Systematic 


6 





4 


8 


6 





2.6 


7.1 


Statistical 
















27.0 


Stat. © syst. 


6 





4 


8 


6 





2.6 


27.9 



Sub-channel 



A u /A s 



WZ — » eve + e~ 
WZ — > /iue + e~ 
WZ — > evji + n~ 
WZ — > fj,ufi + fj,~ 

ZZ -> e + e~vv 
ZZ — >• fi + fj,~vv 



2.242 ± 0.025 459336 
1.495 ± 0.023 419069 
1.704 ± 0.027 493202 
1.443 ± 0.023 443869 

1.638 ± 0.049 319797 
2.052 ± 0.059 342603 



TABLE VIII: Table of acceptance ratios for the different 
sub-channels, where the quoted uncertainties are systematic. 
Also shown are the numbers of observed events at the dilepton 
selection stage, f 58 . 



but still consistent with a standard model prediction of 
3.21 ± 0.19 pb. The ZZ cross section is measured to 
be 1.64 ± 0.46 pb, which is also in agreement with a 
standard model prediction of 1.30 ± 0.10 pb. These are 
the most precise measurements to date of the WZ and 
ZZ cross sections in pp collisions at yfs = 1.96 TeV. 
Correcting for the contribution from 7* and Z/j* in- 
terference and combining with a previous measurement 
in the £ + i~£ + £~ channel yields a ZZ cross section of 



1.44™ pb. 

We thank the staffs at Fermilab and collaborating 
institutions, and acknowledge support from the DOE 
and NSF (USA); CEA and CNRS/IN2P3 (France); 
FASI, Rosatom and RFBR (Russia); CNPq, FAPERJ, 
FAPESP and FUNDUNESP (Brazil); DAE and DST (In- 
dia); Colciencias (Colombia); CONACyT (Mexico); KRF 



16 



Sub-channel ft(xlO" 



WZ -> eve + e~ 
WZ ->■ fiue + e~ 
WZ -t evfi + (i~ 
WZ ->■ pkVji + i>r 

ZZ — > e + e~i^ 



0.70 
0.40 
0.66 
0.61 

0.13 
0.33 



0.20 
0.14 
0.17 
0.16 

0.07 
0.10 
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TABLE IX: Table of 1Z values measured for each of the sub- 
channels, where the uncertainties correspond to statistical and 
systematic components added in quadrature. 
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FIG. 9: Comparison of the measured ZZ and WZ cross sec- 
tions with SM predictions, and with previous measurements 
in leptonic final states. The ZZ cross section measured by 
DO in the ZZ — > £ + £~£ + £~ channel has been corrected to the 
same dilepton invariant mass range as considered in this anal- 
ysis. 



